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“摘要 :在 高 度 计 绝 对 定 标 实验 中 ,利用 差分 GPS(Differential GPS,DGPS) 浮 标 进行 现场 海面 高 测 


( Oit 了 仪器 误差 造成 的 影响 。 前 者 可 通过 低 通 滤波 方式 消除 ,而 后 者 由 于 是 全 频段 的 白 噪声 ,不 能 通 
二 过 低 通 滤波 的 方式 消 除 。 通 过 研究 奇异 谱 分 析 (Singular Spectrum Analysis, SSA) 的 方法 ,确定 了 
,一 谱 入 维 数 和 截断 长 度 的 选取 准则 ,并 对 其 进行 仿真 ,验证 了 该 方法 在 消除 海面 高 序列 的 仪器 误差 吕 
已) 声 中 的 有 效 性 ,同时 也 验证 了 嵌入 维 数 和 截断 长 度 选取 准则 的 有 效 性 。 最 后 ,利用 该 方法 处 理 了 菜 
(次 DGPS 浮标 实验 获得 的 海面 高 序列 ,处 理 结果 与 原始 序列 相 比 较为 平滑 , 且 在 精度 上 有 一 定 的 
加 改善 ,验证 了 奇异 谱 分 析 方法 的 有 效 性 。 
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1s] x 消除 仪器 误差 常用 的 方法 有 最 小 均 方 滤波 和 最 
= 小 二 乘 滤波 四 ,这 两 种 方法 都 要 求 参 考 输 入 序列 与 
-于 在 高 度 计 定 标 中 ,利用 差分 GPS 浮标 获得 的 海 。 原始 序列 相关 。 最 小 均 方 滤波 要 求 输入 序列 必须 是 
面 帘 序列 之 中 通常 包含 了 由 风浪 、 仪 器 误差 等 引起 ”平稳 的 ,而 最 小 二 乘 滤波 是 以 误差 的 平方 和 最 小 作 


的 噪声 分 量 ,需要 对 这 个 原始 时 间 序 列 进行 平滑 处 为 最 佳 准则 的 滤波 方法 ,可 适用 于 非 平稳 序列 。 用 
理 , 从 而 消除 噪声 分 量 , 得 到 星 下 点 的 平均 海面 该 方法 进行 滤波 时 ,需要 知道 序列 中 噪声 的 统计 特 
高 站。 为 了 消除 海面 高 序列 中 的 噪声 分 量 ,通常 的 性 ,从 而 将 一 个 与 原始 序列 中 噪声 的 统计 特性 相同 


做 法 是 对 该 序列 进行 滑动 平均 或 Vondrak 滤波 处 的 序列 作为 参考 输入 ,实现 滤波 。 


理 以 消除 其 高 频 分 量 。 然 而 ,分 析 原 始 序列 频谱 可 由 于 一 个 含有 噪声 的 和 矩阵 的 特征 值 反映 了 有 用 
知 , 由 风浪 引起 的 噪声 其 频谱 仪 集中 在 某 一 频段 内 ， 言 号 和 噪声 的 分 布 情况 ,因此 可 用 奇异 谱 分 析 法 对 序 
而 仪器 误差 引起 的 噪声 分 量 则 分 布 于 整个 频带 上 。 ， 列 进行 去 品 。 利 用 奇异 谱 分 析 法 进行 滤波 时 ,不 需要 
滑动 平均 或 者 Vondrak 滤波 本 质 上 是 一 种 低 通 滤 任何 原始 序列 中 噪声 的 统计 特性 , 仅 需 知道 原始 序列 
波 , 因 此 可 用 于 消除 风浪 噪声 ,但 仪器 噪声 则 不 能 通 大概 的 信 噪 比 即 可 。 而 在 消除 海面 高 序列 仪器 噪声 


过 这 两 种 方法 完全 消除 。 实验 中 ,并 不 知道 任何 仪器 误差 的 先 验 统计 特性 , 


收 稿 日 期 :2014-03-01; 修 订 日 期 :2015-05-24 
基金 项 目 : 中 国 科学 院 国 家 空间 科学 中 心 五 个 重点 培育 方向 (Y32114FB9S) 。 
作者 简介 : 叶 沛 (988 一 ), 男 ,贵州 贵阳 人 ,硕士 研究 生 , 主 要 从 事 高 度 计 定 标 方面 的 研究 。E-mail: yepei1988@163. com. 


H. 


通讯 作者 : 许 可 (1967 一 ), 男 ,辽宁 盖 州 人 ,研究 员 ,主要 从 事 星 载 雷 达 高 度 计 系统 技术 、 合 成 孔径 雷达 高 度 计 系 统 技术 和 信号 处 理 技术 


方面 的 研究 E-mail: xuke(à mirslab. cn, 


662 jun 


仪器 的 绝对 误差 量 与 真实 的 测量 值 相 比 较 小 ,此 时 若 
进行 最 小 二 乘 滤波 ,由 于 没有 任何 仪器 误差 的 先 验 信 
息 ,可 能 因此 带 来 额外 的 误差 。 因 而 可 用 奇异 谱 分 析 
的 方法 对 混 有 仪器 噪声 的 海面 高 序列 进行 降 噪 处理 。 
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2.1 奇异 谱 分 析 基 本 概念 
奇异 谱 分 析 是 一 种 基于 主 成 分 分 析 的 序列 处 理 
Jrik ^ , 它 仅 需 较 少 先 验 信息 即 可 对 原始 序列 进行 
降 品 处理 ,分 析 序 列 结构 (提取 原始 序列 中 的 趋势 、 
周期 成 分 等 ), 因 此 广泛 应 用 于 气象 或 地 理 方面 时 间 
序列 的 处 理 “。 
由 于 和 矩 阵 的 特征 值 能 够 反映 和 矩阵 的 信息 ,并 且 
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可 通过 这 些 特征 值 及 其 对 应 的 特征 向 量 来 恢复 该 矩 
阵 ,大 的 特征 值 对 应 的 是 信号 中 能 量 大 的 组 分 ,小 的 
特征 值 则 可 认为 对 应 的 是 信号 中 的 噪声 信息 。 
此 ,当选 择 用 较 大 的 特征 值 及 其 对 应 的 特征 向 量 来 
恢复 矩阵 时 ,得 到 的 就 是 含有 较 少 噪声 或 不 含 噪声 
的 矩阵 ,这 就 是 奇异 谱 分 析 去 噪 的 基本 思想 。 奇 异 
谱 分 析 分 为 信号 的 分 解 与 重建 两 个 阶段 。 

假设 混 有 了 噪声 的 序列 v; G — 1,2,…,NN) 为 经 
过 了 中 心 化 的 1xX N 的 序列 , 即 序列 x; 的 均值 为 0。 
为 了 能 够 利用 SSA ,可 通过 般 入 的 方法 ,将 序列 uv; 
扩展 为 在 时 间 上 滞后 排列 的 MX (N 一 MM 十 1) 的 空 
[R] A^ s XB RE (trajectory matrix)X。 由 式 (1) 可 以 看 
出 ,矩阵 X 是 Hankel 矩阵。 


Tı T2 Xu TN-M+1 
X2 X3 Xu X N-M+2 

Ac : = [XK Xs Ay nt | (1) 
TM TMH Ti+M TN 


其 中 :M WARA KIAR 
随后 ,计算 空间 坐标 和 矩阵 的 MXRM 自 协 方差 


HERE T, ,这 是 一 个 Toeplitz 矩阵 , 即 任 何平 行 于 主 
对 角 线 上 的 元 素 均 相同 。 


c(0) c(1) c(M — 1)] 
c(l) cO) cCD 
c(l) 
T= : (2) 
c(l) 
IcCM— 1) c(l) cCO) 
其 中 :CG) 是 延 时 为 7 时 xz; 的 协 方差 。 定义 第 个 主 成 分 (Principal. Component, 
| | ON l PC)a; HEREA cG = 1.2. N) TE E E 
G) = = tT ttj 0 x «C M—1 (3) » 
O) = Nm 0«j 的 投影 : 
由 式 (3) 可 知 T. 是 非 负 的 ,对 T. 分 解 特征 值 ， M 
at =X E = rE 1<i<M-1 00 
j=1 


从 大 到 小 的 顺序 重新 排列 ,X41 > A >e 0 Au m0. 
这 些 特征 值 对 应 的 特征 向 量 为 E^, 


( 1 M 
M22 5Fi 
1 M 
AAN — MatE: 
tX 
1 M 
ATE GE; 
aTe Pu 


原始 序列 元 素 x; 下 的 所 有 RC 之 和 即 可 完全 重 
S uis 


RIE at 后 , 即 可 用 来 恢复 原始 序列 ,得 到 重建 
序列 (Reconstruction Component, RC) 9, 


M<i<N-M+1 


ELL M=] (5) 


N—M+2<i<N 


由 于 特征 值 能 够 反映 序列 的 特征 成 分 ,也 就 是 
说 ,利用 特征 值 及 其 对 应 的 特征 向 量 得 到 的 PC 能 
够 反映 序列 的 特征 ,前 面 的 一 些 比较 大 的 PC £63 
序列 中 有 用 信和 号 的 特征 成 分 ,而 后 面 的 一 些 比较 小 
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的 PC 则 主要 反映 噪声 的 影响 。 
此 外 ,对 序列 进行 奇异 谱 分 析 相 当 于 将 序列 按 
其 组 成 成 分 分 解 成 了 若干 份 ,每 个 成 分 的 能 量 都 有 
一 个 特征 值 与 之 对 应 ,并 且 特 征 值 的 大 小 反映 了 各 
成 分 能 量 之 间 的 相对 关系 ,而 序列 的 能 量 是 一 定 的 ， 
Pr R A ER M 越 大 ,特征 值 就 越 小 ,各 个 特征 值 
原始 序列 的 分 解 成 分 就 越 多 ,各 成 分 的 能 量 也 越 小 。 
因此 , 舱 入 维 数 相当 于 分 解 组 成 成 分 的 “分 辨 率 ”。 
利用 SSA 对 原始 序列 进行 去 噪 处 理 时 ,可 以 截 
取 前 面 比较 大 的 特征 值 对 应 的 PC 来 重建 原始 序 
列 , 从 而 达到 去 噪 的 目的 . 设 利 用 前 之 个 PC 恢复 的 
序列 2; 为 : 


p 
Zi = Dr,i=1,2,N (7) 


因此 ,利用 奇异 谱 分 析 的 方法 对 序列 进行 去 品 处 
理 的 关键 在 于 选择 合适 的 庶 和 人 维 数 M 和 截断 长 度 p. 
2.2 骨 入 维 数 及 截取 长 度 的 选取 

设 序列 Yi YI? 分 别 是 对 原始 序列 进行 奇异 
谱 分 析 时 按 其 组 成 成 分 分 解 后 得 到 的 表示 其 两 个 不 
同 成 分 的 序列 .这 两 个 序列 的 加 权 相 关 系数 为 pi 。 
通常 利用 加 权 相 关 系数 oi. 来 评价 两 个 不 同 成 分 的 
序列 间 的 分 离 度 。 


pu^ (YH NT, 
T£ 
| a l Y Ilo 


其 中 : YP |o = VYP YP); (YP ,YD ), = 


N 
N uyy? Gj = 1,2) ju, = min(k,M,N—k +1), 
k=1 


(<4). 


(8) 


2 

式 (8) B ot? 越 小 ,表明 两 个 序列 的 加 权 正 交 性 
越 好 ,它们 间 的 分 离 度 也 越 好 。 若 pi2 = 0, 则 表示 两 
个 序列 完全 加 权 正 交 , 它 们 是 "分离 ” 的 "” 。 

选取 合适 舱 入 维 数 是 为 了 使 奇异 谱 分 解 时 得 到 
的 各 个 独立 的 组 成 成 分 序列 间 的 分 离 度 尽 可 能 高 ， 
以 达到 将 序列 中 的 有 用 信息 完全 从 噪声 中 分 离 出 来 
从 而 对 序列 去 噪 的 目的 .在 DGPS 浮标 得 到 的 海面 
高 序列 中 ,通常 还 包含 了 潮汐 分 量 等 ,这 在 高 度 计 定 
标 实验 中 是 有 用 的 .由 于 藤 入 维 数 相当 于 分 解 得 到 
的 组 成 成 分 的 “分 辨 率 ”, 若 选择 的 嵌入 维 数 过 小 , 则 
测 汐 分 量 等 的 组 成 成 分 有 可 能 与 仪器 噪声 的 组 成 成 
分 部 分 或 完全 混 又 ,使 得 在 对 海面 高 序列 去 噪 时 将 潮 
汐 分 量 等 有 用 信息 也 消除 了 ,导致 定 标 结果 的 准确 度 
和 精度 变 差 ,Hossein"'" 在 数学 上 证 明了 嵌入 维 数 M 


最 大 只 能 取 到 | TS | ([。] 为 ， 的 整数 部 分 ), 且 当 
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Di 
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p 的 选择 和 噪声 信号 与 有 用 信号 之 间 能 量 的 相 
对 大 小 有 关 , 若 p 选 的 合适 , 则 前 个 主 成 分 表示 有 
用 信号 ,而 p 以 后 的 主 成 分 反映 了 原始 序列 中 的 品 
声 . 由 于 在 DGPS 浮标 海面 高 序列 中 ,仪器 误差 对 测 
量 结果 的 影响 较 小 , 若 将 该 海面 高 序列 类 比 为 混 有 
仪器 噪声 的 信号 , 则 仪器 噪声 的 能 量 远 小 于 有 用 信 
号 的 能 量 ,因此 表示 噪声 信号 的 p 以 后 的 特征 值 较 
小 ,由 于 仪器 噪声 是 白 噪声 ,其 能 量 分 布 较 均匀 , 因 
此 表示 噪声 信号 的 p 以 后 的 特征 值 变化 缓慢 . 故 利 
用 特征 值 的 变化 率 来 选择 蕉 断 长 度 p. 

定义 特征 值 变化 率 ; 

r= (Ar —àn)/2sk = 1,2, M—2 (9) 

3C C9) P% r «0.001. HI Ige <— 3 时 ,可 认为 
此 时 户 一 人 一 1 为 截断 长 度 。 


3 ”仿真 及 结果 


3.1 奇异 谱 分 析 法 仿真 结果 

考虑 一 个 由 DGPS 浮标 得 到 的 海面 高 序列 
SSH(D, 

SSH) = SSH, 0) + n(QD) (10) 

其 中 :SSHu(2) 为 真实 海面 高 ,n(7) 为 仪器 误差 。 

假定 真实 的 海面 高 为 SSHu (£) — 104-10 — 4 Xx t, 
t 为 时 间 (SSH GO 单位 为 m, 时 间 t 单 位 为 8) ;仪器 误差 
H o, — 0.1 m, 用 均值 为 0 方差 为 的 高 斯 白 噪声 模 
拟 , 需 要 说 明 的 是 ,这 里 仪器 误差 不 仅 包括 了 浮标 的 仪 
器 误差 ,还 包括 了 DGPS 数据 处 理 中 所 引入 的 误差 。 序 
列 采样 间隔 为 1 Hz, 设 定 采样 时 间 为 长 度 为 1 h, 序 列 
长 度 为 N = 3 600, 则 由 浮标 测 得 的 瞬时 海面 高 原始 序 
列 如 图 1 所 示 。 原 始 序列 的 标准 差 cw = 0.1444 m, 


v^ 


bim 


HAHAY 


0 500 1000 1500 2000 2500 3000 3500 4000 
采样 时 间 /As 
e SSH, « SSH 
1 海面 高 序列 仿真 
Fig.1 Simulated SSH series 
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选择 嵌入 维 数 M = 1 800 ,对 序列 进行 奇异 谱 
分 析 . 图 2 是 序列 分 解 后 得 到 的 特征 值 的 变化 率 曲 
线 .根据 判定 准则 ,& 过 9 时 ,变化 率 > 即 小 于 0.001。 
因此 选择 p = 8 作为 截断 长 度 。 


1 


lgr 
1 
Qo 


图 2 特征 值 变 化 率 


Fig.2 The eigenvalues change rate 


利用 上 述 选 定 的 M RI 再 对 原 序列 进行 奇异 
谱 分 析 , 得 到 去 噪 后 的 重建 海面 高 序列 如 图 3, 它 与 
真实 海面 高 序列 SSH, 基本 重合 .重建 的 序列 的 标 
ME o, = 0.1058 m, 


10.15 AP 
10.10 M 
1005[ — d 

10.00 


0 500 1000 1500 2000 2500 3000 3500 4000 
采样 时 间 /s 
，SSH。 : 重建 海面 
3 ”重建 的 海面 高 序列 
Fig.3 The reconstructed SSH series 


应 用 p 以 后 的 主 成 分 来 重建 噪声 信号 ,图 4 是 
重建 的 噪声 序列 的 频谱 ,可 以 看 出 ,其 功率 在 整个 频 
段 上 分 布 均匀 ,符合 白 噪 声 的 频谱 分 布 特性 .图 5 是 
重建 噪声 的 幅度 分 布 直方 图 ,其 幅度 符合 高 斯 分 布 。 
重建 仪器 噪声 的 标准 差 o,, = 0.09958 m, 均 值 m = 
— 5.613 X 10”m, 基 本 符合 添加 的 噪声 的 统计 特性 。 
因此 ,从 仿真 中 可 以 看 出 ,用 奇异 谱 分 析 的 方法 
能 够 较 好 地 将 序列 中 的 有 用 信号 和 噪声 信和 号 分 离 出 
来 ,其 去 噪 效果 较 好 。 
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0.012 


0.010 


0.008 


0.006 | 
| | | | | | H 
0.002 f "E m Is VEIT 

0 0.1 0.2 0.3 0.4 0.5 


频率 /Hz 
图 4 重建 的 噪声 序列 频谱 


Fig.4 The frequency spectrum of reconstructed noise series 
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出 现 次 数 


-04 -03 -02 -01 0 01 02 03 04 
仪器 误差 /m 
图 5 仪器 噪声 幅度 分 布 直方 图 


Fig.5 The distribution histogram of instrumental noiseamplitude 


3.2 WIE ERAN 2E A K REGE DS ER s R AS E TI 
为 了 进一步 比较 不 同 的 M 和 p 对 数据 处 理 结果 
的 影响 ,同时 为 了 与 上 节 仿 真 结 果 比 较 , 分 别 用 固定 
p 二 8 改变 M 和 固定 M== 1 800 改 变 p 这 两 种 方法 处 
理 上 述 海面 高 序列 。 表 1 是 p 二 8 时 不 同 M 下 得 到 的 
字 列 的 标准 差 , 表 2 是 M = 1800 时 不 同 p 下 得 到 的 
序列 的 标准 差 。 由 奇异 谱 分 析 去 噪 的 原理 可 知 ,p 必 
须 小 于 M.H p 宇 1, 故 在 表 1 中 MM 从 M = 9 变化 到 
M = 1 800; 表 2 中 p 从 p= 二 1 变化 到 p= 二 1 799, 
比较 表 1 及 表 2 可 以 得 到 ,重建 海面 高 序列 的 
标准 差 是 随 着 M 的 增 大 而 减 小 , 随 着 p 的 增 大 而 增 
大 的 ;重建 噪声 序列 的 标准 差 是 随 着 M 的 增 大 而 越 
接近 仿真 中 设 定 的 噪声 标准 差 , 随 着 p 的 增 大 而 越 
远离 设 定 的 噪声 标准 差 的 ,虽然 表 2 中 p= 二 2 时 重建 
海面 高 序列 的 标准 差 最 小 ,但 其 对 应 的 重建 噪声 序 
列 标准 差 较 p = 二 8 时 偏离 较 多 , 故 在 M = 1 800 下 ， 
p= 二 2 不 是 最 优 的 .因此 ,由 以 上 分 析 可 知 ,在 上 节 仿 
真 中 选择 M = 1800, = 8 是 较 合 适 的 。 


第 4 期 "pO 沛 等 :基于 奇异 谱 分 析 的 DGPS 浮标 海面 高 测量 误差 研究 665 
R1 p=8 时 不 同 M 下 得 到 的 序列 标准 差 
Table 1 The series std of Various M with p—8 
M 9 500 1 000 1 500 1 800 
重建 海面 高 序列 标准 差 0. 13566 0. 10696 0. 10657 0. 10635 0. 10581 
重建 噪声 序列 标准 差 0. 03096 0. 09833 0. 09919 0. 09949 0. 09958 
R2 M=1 800 时 不 同 p 下 得 到 的 序列 标准 差 
Table 2 The series std of Various p with M—1 800 
p 1 8 500 1 000 1 500 1 799 
重建 海面 高 序列 标准 差 0. 09490 0. 10583 0. 12874 0. 13914 0. 14444 0. 14562 
重建 噪声 序列 标准 差 0. 10132 0. 09955 0. 06202 0. 03727 0. 01603 0. 00284 
AN NI m «y HH uu 一 
4 利用 SSA 消除 海面 高 序列 仪器 噪声 


下 面 以 某 次 水 库 DGPS 浮标 实验 为 例 , 对 其 进行 
奇异 谱 分 析 来 消除 序列 中 仪器 误差 的 影响 。 由 于 实 
验 当 天 水 面 无 风 或 风速 很 小 ,因此 ,水面 无 风浪 或 者 
风浪 很 小 ,可 认为 浮标 得 到 的 海面 高 仅 受到 仪器 误差 
的 影响 。 查 看 浮标 采集 到 的 GPS 数据 并 截取 其 中 质 
量 较 好 的 一 段 ,利用 GAMIT/GLOBK 软件 中 的 动态 
处 理 模 块 TRACK 解 算 , 得 到 浮标 所 在 位 置 的 海面 高 
数据 。 图 6 是 浮标 在 大 地 坐标 下 高 度 ( 海 面 高 ) 方 向 
的 运动 轨迹 ,其 标准 差 o, —2. 88 cm, 
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采样 时 间 /s 
图 6 浮标 在 HH 方向 运动 轨迹 


Fig.6 The buoy H orientation movement path 


选择 M—1 800 对 序列 进行 奇异 谱 分 析 。 图 7 
是 分 解 得 到 的 特征 值 的 变化 率 。 根 据 判定 准则 选择 
5710, 

利用 上 述 选 定 的 参数 ,再 次 对 原 序 列 进行 奇异 
谱 分 析 处 理 , 得 到 的 重建 序列 如 图 8 所 示 , 与 原 序 
列 相 比 , 重 建 的 海面 高 序列 较 平滑 。 其 标准 差 o 5 
2.69 cm. E RSS 意义 上 ,标准 差 较 人 处理 前 改善 了 
约 1. 03 em, 


0 200 400 000 3800 1000 1200 1400 1600 1800 
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图 7 特征 值 变 化 率 


Fig.7 The eigenvalues change rate 
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”原始 序列 。 。 重 建 序列 
8 重建 的 海面 高 序列 
Fig. 8 The reconstructed SSH series 
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在 高 度 计 绝 对 定 标 实验 中 ,利用 DGPS 浮标 进行 
现场 海面 高 测量 得 到 的 海面 高 序列 中 ,通常 包括 了 仪 
右 误 差 造成 的 影响 。 由 于 在 浮标 布 放 期 间 , 仪 器 误差 
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是 白 噪声 ,因此 不 能 通过 低 通 滤波 的 方式 完全 消除 仪 
器 误差 对 海面 高 序列 的 影响 。 奇 异 谱 分 析 是 一 种 基 
于 主 成 分 分 析 的 序列 处 理 方法 , 仅 需 较 少 的 先 验 信息 
即 可 对 序列 进行 降 品 处理 。 本 文 研 究 了 奇异 谱 分 析 
的 基本 理论 及 如 何 选择 合适 的 府 人 维 数 和 截断 长 度 ， 
并 通过 仿真 验证 了 奇异 谱 分 析 方 法 及 骨 人 维 数 和 截 
断 长 度 选取 准则 的 有 效 性 。 随 后 ,利用 奇异 谱 分 析 处 
理 了 某 次 水 库 DGPS 浮标 实验 的 海面 高 序列 ,处 理 结 
果 在 精度 上 有 一 定 的 改善 。 因 此 奇异 谱 分 析 法 在 消 
除 海面 高 序列 仪器 噪声 中 较为 适用 。 
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Study on Denoising the Instrumental Errors of the Sea Surface 


Height Series Derived from Buoys 
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Abstract: In absolute altimeter calibration campaigns. the Sea Surface Height (SSH) series derived from the in-situ 


Differential GPS (DGPS) buoys always include high frequency errors such as wind waves and instrumental noises. 


The former can be removed by low pass filters , while the latter can not, because the noises are white and can spread 


in the whole frequency band of interest. In this paper. the Singular Spectrum Analysis (SSA) method was applied to 


solve the problem,and select criteria of the embedding dimension and intercepting length. Simulations based on this 


innovative method were executed, the effectiveness in removing the instrumental noise and the selction criteria was 


validated. Finally, SSH series derived from an in-situ DGPS buoy experiment were processed by this method, and 


the study suggested that the processed series were smoother compared to the original series, and had some improve- 


ments in accuracy. Therefore, the validity of the SSA method was demonstrated. 


Key words: Altimetric calibration; SSH series; Instrumental errors; Singular Spectrum Analysis(SSA) ; De-noise 


